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STABILIZING  EFFECT  OF  GAS  CONDUCTIVITY  (VOLUTION  ON  THE 
RB8IOTiyE  SAUSAGE  MODE  OF  A  PROPAGATING  BEAM 


1.  Introduction 

A  self-pinched  oloctron  or  ion  bona  propagating  in  gas  will  excite  a 
substantial  return  currant  if  tha  gas  is  prs-ionissd,  if  it  is  rapidly,  ionised 
colllsloaally  by  tha  baaa  hsad1>  or  if  avalanche  breakdown  is  driven  by 
Inductive  electric  fields  at  the  beau  head*  One  nay  well  expect  highly 
current  neutralised  beans  to  be  subject  to  a  variety  of  instabilities  excited 
by  the  repulsive  force  between  bean  and  return  current*  and  several  nodal 
calculations^*^  have  reached  this  conclusion  with  regard  to  beans  propagating 
in  resistive  plasna*  The  situation  differs  sharply  fron  that  of  a  non- 
current-neutralised  bean  in  a  resistive  plasna,  for  which  the  hose  node  is  the 
only  unstable  node* 

Under  conditions  where  they  are  unstable,  the  axisynnstrlc  bean  nodes 
appear  to  be  particularly  dangerous  to  propagation,  because  they  are  alnost 
inevitably  excited  at  large  anplltude.  Unless  the  bean  enittance  is  perfectly 
Hatched  at  injection,  the  baaa  will  oscillate  in  radius,  and  in  particular, 
the  violent  plnchdown  associated  with  the  process  of  nose  expansion  and 
erosion* »•  can  be  expected  to  excite  sons  radial  oscillations*  Non- 
axlsyanetrlc  nodes,  such  as  hose  and  fllaaantatlon,  nust  grow  out  of  initially 
low-level  noise  if  the  accelerator  produces  a  high-quality  bean,  and  thus  nust 
e-fold  aany  aore  tines  before  they  pose  a  threat  to  bean  Integrity. 

We  classify  the  linearised  nornal  nodes  of  a  bean  by  an  asinuthal  node 
nunber  n  and  a  radial  node  nunber  a.  For  a  node(a,n),  all  perturbed 
quentltlea  ♦  are  of  the  fora 

f(r,9,s,0  -  exp(la0  -  IQs)  J(r,0, 

where  c  •  ct  -  s  is  the  poeition  in  the  been  neesured  beck  fron  the  been 

A 

heed,  loughly  specking,  n  is  the  nunber  of  oscillations  of  t(r,()  as  r  varies 


froa  sero  to  th«  outside  of  the  beaa,  for  fixed  C  •  The  a  ■  0  (axisyaaitrlc) 
•odes  ere  soaetlaes  referred  to  collectively  as  sausage  nodes ,  but  we  shell 
reserve  this  tera  for  the  a  ■  0,  n  -  1  node,  which  corresponds  roughly  to 
self-siallar  radial  oscillation  or  "breathing"  of  the  beaa,  and  we  shall  refer 
to  the  next  axiayaaetric  aode,  a  -  0,  n  -  2,  as  the  axlsyamtric  hollowing 

4 

■ode.  A  recent  theoretical  study  by  Oha  and  Laape  ,  based  on  the  slapllfying 

assuaptlons  of  a  fixed  flat-topped  radial  profile  of  plasaa  conductivity  and  a 

flat-topped  beaa  radial  profile,  predicted  sausage  inetablllty  when  the  ratio 

of  plasaa  return  current  to  beaa  current ^  I  /I.  >  0.50  and  axisyaastrlc 

P  6 

hollowing  instability  when  1  /l.  >  0.38.  Lee3  extended  the  theory  of  the 

P  b 

sausage  aode  to  arbitrary  beaa  profiles  and  arbitrary  (but  fixed)  conductivity 

profiles,  and  also  predicted  sausage  instability  when  I^/I^  exceeds  a 

threshold  whose  exact  value  depends  on  the  profile.  For  similar  Bennett 

profiles  of  beaa  current  and  conductivity,  the  threshold  is  0.69;  when 

I  / I  >  0.75,  instability  at  fl  •  0  is  predicted.  However  the  sausage 
P  b 

instability  has  not  been  observed,  as  far  as  we  are  aware,  for  beaas 
propagating  in  neutral  gas,  although  I^/I^  often  exceeds  any  of  these 
instability  thresholds. 

In  this  paper,  we  present  a  aore  coaplete  linearised  theory  of  the 
sausage  aode  of  a  relativistic  electron  beaa  which  begins  froa  Lee's 
foraulation*  but  Includes  a  self-consistent  treataant  of  the  plasaa 
conductivity,  including  collislonal  ionisation  of  tha  gas  by  the  perturbed 
beaa.  We  find  that  the  conductivity  channel  perturbs  in  such  a  way  as  to 
follow  the  distortions  of  the  beaa.  This  inhibits  the  separation  of  beaa  end 
plasaa  current,  which  is  the  cause  of  the  instability,  and  consequently  leads 
to  a  aore  stringent  lnsteblllty  condition,  which  is  never  satisfied  for  besas 
propagating  in  high  density  initially  unionised  gas.  (The  situation  at  low 


iu  density,  <  50  torr  depending  on  conditions,  Is  discussed  Inter  In  this 
introduction.)  This  instability  condition  Is  [see  Eq.  (41)  In  See*  2) 

(i> 

•ft .  A 

where  is  s  radially-averaged  net  current  which  determines  the  assn  pinch 
force. 
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Klb  ,  d  /*o( r-0,c)s2x 

tt  -  ar ( — is - ' 


is  s  esssure  of  the  rate  of  change  of  conductivity  a  due  to  beam-molecule 
collisions,  k  Is  i  coefficient  depending  on  the  ionisation  coefficient  and 

mobility  of  the  particular  gas,  and  H  Is  a  factor  of  order  unity. 

> 

The  inequality  (1)  can  easily  be  satisfied  for  beaas  propagating  Into 
pre-lonlsed  gas.  The  Instability  predicted  under  those  conditions  typically 
has  a  smaller  growth  rate  than  the  hose  node,  but  could  dominate  if  the 
sausage  node  is  initiated  at  larger  amplitude,  as  it  normally  will  be  in  a 
well-prepared  beam. 

For  beans  Injected  into  neutral  gas,  a  very  brief  burst  of  avalanche 
ionisation  at  the  beam  head  typically  has  a  strong  Influence  on  the  degree  of 
current  neutralisation*  This  effect  is  tacitly  included  in  our  theory,  since 
lagf/lb  is  treated  as  a  free  parameter.  However  our  analytic  theory  treats 
only  beasrcolllslonal  ionisation,  and  not  avalanche,  in  the  beam  body  Hhere 
the  instability  grows.  This  is  usually  a  good  approximation,  except  in  low 
density  gas  (5-50  torr  depending  on  beam  current  density  and  gas  type).  He 
have  found  some  cases  in  low  density  gas  where  avalanche  at  the  beam  heed  is 


so  strong  that  the  instability  condition  (1)  is  satisfied,  but  la  all  of  these 


caees  noted  to  date  avalanche  should  also  ba  locludad  in  modeling  in  tha  haaa 
body.  Heuristic  eonsidarationa  indicata  that  tha  naglactad  avalanchs  tar*  in 
tha  baa*  body  would  snhanca  stability,  by  furthar  parturbing  the  conductivity 
channal  so  as  to  follow  tha  baa*  distortions,  We  have  tasted  this  idea  by 
performing  a  few  simulations  with  tha  axlsyaaatrlc  baa*  envelope  code  VIPER-0, 
which  includes  avalanche  everywhere  and  peraits  sausage-like  oscillations  to 
develop  (but  does  not  perait  any  higher  axlay*aetric  perturbations).  Sausage 
instability  did  not  occur  in  these  cases,  even  though  condition  (1)  was 
satisfied.  Further  study  of  this  low-density  regies  is  needed,  however. 

Particle  siaulations  at  several  laboratories8-10  have  recently  observed 
strong  axisy**etrlc  Instabilities  under  a  variety  of  conditions  that  are 
coapatible  with  the  instability  condition 

V1*  *  <2> 

but  are  incompatible  with  condition  (1).  We  have  shown  by  naans  of  a 
slaula t ion  analysis  that  these  instabilities  Involve  the  hollowing  node,  not 
the  sausage  node,  and  are  triggered  by  a  couples  set  of  circumstances  with 
other  requirements  in  addition  to  (2).  These  results  are  reported  in  a 
separate  paper.11 

The  outline  of  this  paper  is  as  follows.  We  introduce  our  model  and  list 
ite  aaeuaptions  in  Sec.  2.  Our  analytic  calculation  of  the  sauaage  node 
dispersion  relation  is  presented  in  Sec.  3.  Our  conclusions,  with  regard  to 
beam  propagation  in  initially  neutral  gas,  pre-ioniaed  gas,  or  in  a  channel  of 
fixed  conductivity  profile,  are  discussed  in  Sec.  4. 


2.  Formalism  and  AwmptloM 

la  this  Motion  we  develop  a  fully  a— lytic  theory  of  tha  sausage  mode  of 

a  relativistic  alaetroa  beau  that  Includes  tha  modifications  of  tha  chan— 1 

conductivity  that  result  froa  beae-colllsional  ionization  treated  — If- 

consistently  with  the  sa— aging  of  tha  bean*  He  find  that  the  conductivity 

chan— 1  tends  to  follow  the  —usage  distortion  of  the  beau;  as  a  result* 

spatial  separation  of  the  plas—  return  current  de— ity  J  from  the  beau 

P 

current  de— ity  is  reduced,  and  tha  node  is  found  to  be  such  —re  stable 

than  it  would  be  in  a  fi— d  conductivity  chan— 1.^ 

In  order  to  carry  out  tha  a— lysis  in  siuple  form,  we  —he  a  nuaber  of 
siupllfylng  assumptions,  —at  of  then  having  wide  validity*  Four  of  thea e 
assumptions  have  been  widely  employed  in  beam  propagation  theory*  They  are 
that  the  background  gas  can  be  regarded  as  an  in— bile  medium  with  a  scalar 
conductivity  o(x,t),  that  the  beam  is  highly  relativistic 
(Y  »  1)  and  paraxial  (v  «  v  for  all  electrons),  and  that  therefore 

A  Z 

V  ■  c. 
z 

We  consider  only  1— tablllty  growth  in  the  region  of  the  beam  which  we 
shall  call  the  be—  "body",  behind  the  pi— h  point1’6  but  forward  of  the  beam 
tall  where  recombl— tion  limits  the  conductivity.  This  is  the  region  where 
violent  axisymme trie  Instability  has  been  observed  in  simulations,8”10  and 
where  theory  indicates  that  1— tabllltles  should  grow  — st  rapidly.  Here  the 
conductivity  a  is  large  enough  to  1— ure  space  charge  — utrality  in  the 
vic.i-ity  of  the  be—,  and  Maxwell's  eq— tlo—  reduce  to  Ampere's  law, 

TTi'h  ~  — frfca  *  -  r-  V*-*,.)  <3> 


for  an  axisym— trie  beam,  where  A(r,(,z)  is  the  axial  compo— nt  of  the  vector 


potential.  W«  can  aspect  the  bean  body  in  equilibrium  to  have  a  Bennett 
radial  profile, 


Jbo(r’C»x) 


bo 


,,  2,  2,2 

(1  +  r  /ao) 


(4) 


(Other  profiles  could  be  treated.)  In  the  bean  body,  the  Bennett  radius 

a  typically  decreases  slowly  with  C  and  increases  very  slowly  (erosion  and 
o 

Horde lack  expansion1' »®)  with  s;  we  use  the  approximation  that  a  is  completely 

o 

independent  of  both  s  and  C  in  equillbriua.  Since  avalanche  is  normally 
unimportant  in  the  bean  body  except  for  beams  propagating  in  low  density  gas, 
e.g.  P  <  50  torr  for  a  typical  induction  linac  beam  with  1^  ■  10  kA, 
a^  •  0.5  cm,  and  since  recombination  is  by  definition  unimportant  in  the  beam 
body,  we  assume  that  the  conductivity  evolution  is  determined  by  direct 
colllslonal  ionisation  of  the  gas  by  the  bean. 


Sc 


o(r,c,s)  -  kJ, 


b’ 


(5a) 


(However  the  effect  of  avalanche  in  the  beam  head  is  tacitly  Included  as  an 

initial  condition,  l.e.  the  values  of  o  and  Ip/I^  specified  at  the  front  of 

the  beam  body  segment  Included  in  the  calculation,  C  •  C  ,  depend  implicitly 

o 

on  avalanche  in  the  beam  head,  C  <  CQ»)  We  neglect  the  fairly  weak  dependence 
of  o  on  plasma  temperature  (through  the  plasma  electron  collision  frequency  u^) 

end  thus  treat  k  as  a  constant.  A  reasonable  estimate  of  k,  based  on  the  typical 

12  12  -1 

value  v  *  2  x  10  sec  for  air  at  standard  density  and  0.5  eV  <  T  <  1  eV,  is 

Q  © 

k  •  8,8  x  10’U  cm/statcoul.  (5b) 
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To  facilitate  th«  analysis,  vs  also  asks  ths  following  additional 

asaumptions.  Ths  equilibrium  electric  field  E  (c,s),  which  has  a  weak  radial 

so 


dependence 


2  2 
1  +  b4/a‘ 

E*o(r»5»z)  *  tn  C - yy)» 

1  +  r  /a* 
o 


Is  regarded  as  Independent  of  x, 


Eso(r,5*E)  “  ^o^0’5***' 


In  (6a)  b  Is  a  large  radius  where  charge  neutrality  falls.  Approximation  (6b) 
Is  common  to  all  previous  stability  analyses.  He  note  elsewhere  that  the 
breakdown  of  this  approximation  plays  a  key  role  in  the  destabilisation  of  the 
hollowing  mode, I*  but  we  believe  that  the  approximation  is  acceptable  in 
analysis  of  the  sausage  mode. 

As  a  result  of  Eqs.  (4)  -  (6),  the  equilibrium  conductivity  and  plasma 

current  density,  o  (r,0  and  J (r,C)»  both  have  Bennett  profiles  of  radius 

o  po 

a0«  This  turns  out  to  be  very  helpful  to  the  analysis:  it  allows  us  to 
reduce  a  problem  in  r,C  and  s  to  the  form  of  ordinary  differential  equations 
in  c  and  s  only.  For  related  problems  of  interest,  the  independence  cannot  be 
eliminated.  For  example,  in  the  beam  "tail”  where  beamrcolllslonal  ionisation 
is  balanced  by  recombination,  o(r)  «  [ J^(r)]  ,  a  broader  profile  than 

J^(r).  As  a  result,  much  of  the  plasma  current  flows  outside  the  beam  and  has 
no  destabilising  effect  on  the  beam,  Which  essentially  guarantees  that  the 
sausage  mode  will  be  stable.  On  the  other  hand.  If  the  channel  Is  fully 
Ionised,  as  is  the  case  typically  for  applications  to  Ion-beam  inertial 
fusion,  o(r)  ■  t  3/2  m  [t  (r)]3/2  for  spitser  conductivity  and  no  effective 


.  w  4-  *-■*-* * 


heat  loss  mechanism.  Thus  o(r)  iw  narrower  than  J^(r) ,  causing  the  plasma 

current  to  peak  on  axis,  further  destabilizing  all  bean  nodes. 

We  assune  for  convenience  that  and  Jpo  are  also  Independent  of 

(The  analysis  could  be  carried  out  without  these  assumptions,  but  not  in 

closed  forn.)  The  Justification  for  the  latter  is  that  the  decay  length  t  of 

o 

Jfeo  1*  long  conpared  to  the  decay  length  of  current  perturbations, 


t  .  *231  t« 
o  c  a 

o 


1  2c 


also  characterizes  the  instability  growth  length. 

tfe  consider  only  beam  perturbations  in  the  form  of  self-similar  radial 
expansion  and  contraction,  i.e.  the  perturbed  bean  is  of  the  form 


Jb(r,C,z)  -  — j 


wa  (C,z)  [l  +  r  /a  (C,z)]‘ 


which  is  a  reasonable  but  not  exact  representation  of  the  sausage  mode.  The 
variation  of  the  root-nean-square  bean  radius*^  a  (€,z)  is  calculated  fron  the 
Lee-Cooper*-*  envelope  equation, 


a2;  ,  c2  _  0^ 

"at7  a7  a 


where  U*  is  a  measure  of  the  average  pinch  force. 


,  21  2wrJ.  (r)  ,  2*r  [j.  (r  )  +  J  (r  )] 

U2h  -=*  /  dr  — =-£ -  /  rdr  - ^ - 2 - - 

lA  o  ib  o  S 
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I  ■  17f  kA  is  Che  Alfven  current,  J  is  the  plasma  current  density^,  and  e  is 

A  r 

the  emlttance ,  defined  as  a  <9>,  where  <8>  is  the  roo  t-mean-s quare  beam 
electron  velocity  angle. 

The  principal  effect  omitted  by  the  model  (8)  -  (10)  is  sausage 

oscillation  damping  due  to  phase  mixing  among  beam  electrons  of  different 

I 

betatron  frequency.  This  effect  is  included  phenomenologically  by 

adding  a  damping  term  in  the  form  derived  by  Lee  and  Yu^* 


l£_ 

dz 


(ID 


The  damping  constant  a  is  sensitive  to  the  beam  profile. ^  Lee^*^  estimates 
a  -  0.7  for  use  with  a  Bennett  profile. 

We  perform  a  linearized  perturbation  calculation  in  which  small 


perturbations  are  added  to  the  equilibrium  profiles  of  Jb,  Jp» 
o  and  Ag  described  above.  Perturbed  quantities  are  calculated  from  the 
linearized  forms  of  Eqs.  (3),  (5),  (9),  (11),  and  Ohm's  law, 


JmE  .  (12) 

P  * 

Because  of  the  assumed  radial  dependences  of  the  equilibrium  and 
perturbations,  the  equations  reduce  to  coupled  ordinary  differential 
equations. 

We  note  particularly  that  both  the  equilibrium  and  the  perturbed 


conductivity  are  calculated  self-conslstently  from  Eq.  (5);  this  is  the  new 
feature  of  the  present  calculation  and  leads  to  enhanced  mode  stability* 


3.  Calculation 


The  equilibrium  described  in  Sec.  2  is  specified  by' 


Jbo-  \o(1  +  rX>'2> 


(13.) 


~  2.2-2 
J  -  J  (1  +  rVaZ)  , 
po  po  o 


(13b) 


o  -  kX  C(1  +  r2/a2)-2, 

O  DO  O 


(13c) 


1  +  r2/a2  _  2 

V  v°*n - 2TT  "  ~  V°  in  T  • 

1  +  b  /a*  a*" 

o  o 

dA  dA  w2 
_  _ o  o  „  b 

Eor“  “  d?  dC  a2  * 


(13d) 


(13e) 


Equations  (3),  (9)  and  (10)  reduce  to  the  equilibrium  relations 


dA  ira 


A  +  t  -j-r 
o  o  dC 


-  7^  Jbo* 


(14) 


2.  -2  2  _ 

c  /  a  -  U  •  I  C,/I., 
o  o  o  eff  A’ 


(15) 


where 


eff  5 


(16) 


for  the  present  case  of  beam  and  plasma  currents  with  the  same  profile2 

The  front  of  the  beam  segment  of  interest  Is  taken  to  be  at  C  ■  C  ;  the  value 

o 

of  may  be  chosen  so  that  Eq.  (13c)  gives  the  conductivity  desired  as  an 
initial  condition  in 

10 


The  perturbed  pert*  of  Jp,  o  and  A,  designated 
SJb,  5jp,  do  and  5A,  ere  treated  as  snail  quanltles.  To  first  order  Eqs.  (3) 
(5),  (9),  (11)  and  (12)  bacons 


13  ^  3  tk  4™o  36A  4ir  3Ao 
r  Tv  r  77  6A  "  c”  7T  c“  TT 


JSi  -  SSL  ®  do  -  SS,  6jt 
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(17a) 
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where  we  have  used  (11)  and  (15)  to  simplify  (9).  These  equations  adnit  the 
solution 


i  -  A.* 

5J."  ”  Jv((|l)  j  j  i  » 

b  b  (l  +  r2/a2)3 


(19a) 


"  J  U,*) 


1  -  r2/.2 
_ o 


P  (1  +  r/.y 


(19b) 


-  -  o(C, *) 


-  A  (5,*) 


,  2,  2 
1  ~  v*t , 

(1  +  r2/.2)' 

)i^£ 

1  +  r  /a 


(19c) 


(I9d) 


which  corresponds  to  self-slnllar  oscillations,  since 
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(20) 


1  -  r2/.J 

- 2 — FT 

(1  +  r2/a2)3 
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Equation*  (17)  and  (18)  than  raduca  to  tha  ordinary  diffarantial  aquations 


2  ~  2  2  ~ 
r,  *V  Jbo  .  3  'i  ?  _  wao  ‘■j  .  »*o  “o  fb  i  : 

(1  + - j- - ?  55 J  A  “  V  “2?  dT  tn  °* 


(21a) 
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(21d) 


We  address  the  linearized  envelope  equation  (2ld)  first.  Since  1^  and  Ip 

are  unchanged  by  aelf-sinilsr  expansion,  only  the  cross  terns  between  Jp  and 

2 

Jb  in  Eq.  (10)  contribute  to  the  perturbed  pinch  force  60  .  Recalling  our 
convention  that  Ip  and  1^  are  both  positive^,  we  find 


6U2  _  1  ,  \  ~  ,  h  x 

~2~  "  t-  i~  Jb+  I - V  » 

0  3J.  eff  eff  p 

O  DO 


(22) 


and  the  perturbed  envelope  equation  (21d)  then  reduces  to 


323.  33.  ala  xk 

*  --„j>  -(2-4_e)i+4  ‘  j. 

)  3  x«ff  b  3  x-#f  P 


3(*Ag) 


TOT 
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eff 
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where  we  have  defined  an  average  betatron  wavelength 


V  W 


(24) 
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-v 


The  quantity13  a  appears  only  through  X  . 

o  p 

In  order  to  eloaa  the  analyela  t  we  auat  express  3  in  teraa  of  3.  though 

p  a 

Sqe.  (2U-c).  first  we  simplify  Bq.  (21a)  by  using  the  approximation  that  Ip 
as  well  as  1^  are  independent  of  c,  eo  that7 


-  *  const 


and  from  Eqe.  (12)  and  (13e) 


dA  .2 

Jpo"  K  3boc  dTla7  * 

*o 


y  s 

k  x 


(The  analysis  could  be  carried  out  without  these  approximations ,  but  would 
then  require  numerical  solution  of  a  complicated  ordinary  differential 
equation,  rather  than  yielding  the  solutions  we  shell  find  in  dosed  fora*) 
Thus  Bq.  (21a)  can  be  written  as 


(IfL  +  * 


I  ~ 


1  r>  11.1  _  P  O 

JboC  7Z}  b  r 


Next  we  rewrite  Bq.  (21b)  as 


V  (M) 

and  use  (28)  to  ellaate  o  from  (27),  which  yields  an  expression  for  A  in  tei 
of  Jb, 

U  ♦  c  *  *  3boc  1?^  *  "  +  c  TT  "  if)  "V  <2*> 


We  also  use  Bqs*  (13),  (19)  and  (28)  la  (17e),  to  obtain 


V m  t  "  *  Jw  '  Tc  •  <J*) 

and  ulat  (28)  again  la  (30)  giant 

<3l> 

Equation*  (29)  and  (31)  togathar  specify  Jp  a*  a  llnaar  function  of  Jfe. 

Ua  obaarva  that  (23),  (29)  and  (31)  conatltuta  a  closad  systea  of  llnaar 

aquations  with  constant  coafflciants  if  (3/3s)  and  C(3/3c)  =  (3/9  In  C)  ar* 

ragardad  as  tha  basic  operators.  Thus  It  is  natural  to  usa  s  and  In  C  as  tha 

Indapandant  variables.  A  couplet*  solution  of  tha  problna  with  perturbation 

initial  conditions  at  c  a  C  would  require  a  Laplace  transfora  analysis,  which 

o 

is  beyond  tha  scope  of  this  paper,  but  tha  aquations  adult  a  free  node 
solution  In  tha  fora 


f(C,*)  -  f  exp  [-ifla/Xg-  itf  in  (C/CQ)] 


„  -lfls/X 

•  fa 
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(32) 


where  f(C,s)  is  any  perturbed  quantity.  Tha  * -dependence  of  the  node  is  In 

tha  usual  exponent ial  fora,  but  in  ?  tha  aoda  shows  power  law  growth  or  decay 

and  oeclllatloa  with  steadily  Increasing  wavelength,  due  to  tha  non-uniforalty 

of  tha  equillbrlua  in  (,  l.a.  tha  llnaar  increase  of  e^(r,5)  with  {. 

Using  (32)  In  (29)  and  (31)  yields  tha  required  relation 

between  3  and  3  , 
p  o 


3  -  [- 

*  i  -  if 


UU  -  -  ytt) 

x  i  Jb* 
r  -  la  ) 


(1  -  la)( 


TT 


-  •: 


(M> 


1A 


We  also  dm  (32)  to  reduce  the  envelope  equation  (23)  to  algebraic  fora,  and 


uaa  (33)  to  eliminate  J  .  After  eoae  algebra,  but  no  approximations,  the 

P 

envelope  equation  reduces  to  the  diaperelon  relation 


-  0  -  iafl  -  -  2 
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where 


X 


-  fli-Ii 

'  d;  "2c 
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(35a) 


from  Eq.  (7b)  Is  the  characteristic  decay  tiae  of  the  perturbed  plasaa 
current,  and  o  le  the  phenoaenologlcal  damping  constant  froa  Eq.  (11). 

If  we  use  Eq.  (5b)  for  k,  we  can  write  Eq.  (35a)  in  the  convenient  fora 


4.  Iwulu 

The  dispersion  rslscion  (34)  can  ba  solved  for  either  Q  or  m  as  a 
function  of  the  other.  The  node  is  unstable  if  Q^>  0  for  real  w 
or  if  «i>  0  for  real  Q.  We  consider  the  instability  condition  Q^>  0  for  real 
Z.  Re writing  (34)  in  the  schanatic  form 


-  a*6  -  loft  -  f^x,  ip/i,££.«)  +  ip/i.££,  S), 


where  and  F2  are  real,  this  condition  is 


Ft  >  -  F22/o2. 


For  a  *  0  (no  phase  nix  deeping  included  in  the  foraallsa)  all  nodes  are 
unstable;  as  eight  be  expected  the  conditions  for  instability  bacons  eore 
restrictive  as  a  increases,  but  Instability  can  occur  even  in  the  Unit 
a  ♦  •  if  Fj>  0,  l.e. 


^2 

1  +  u  eff  « 


There  ere  nodes  Z  that  satisfy  condition  (38)  if  and  only  if 


and  the  range  of  unstable  nodes  is  given  by 


Ct-  C2  <  S2  <  Cf  c2. 


(40a) 
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where 


V  “T  Ul  *  x>  t*— “  4x2  *  31» 

1  8X*  off 

c  i-L  {[(u  x>  I* —  4X2  -  j]2  -  ux2}1/2. 


(40b) 


(40c) 


If  condition  (39)  ia  satisfied,  there  aro  vary  firn  nodal-independent 
grounds  for  expecting  strong  instability*  independent  of  tha  il saying 
coafficiant  a.  In  fact*  if  ona  solves  for  «(S)  it  ia  soon  that  instability 
can  occur  (u^>  0)  in  this  case  area  with  0  •  0,  l.e.  asra  tlas-lndepeodeat 


baaa  non-uni f orui t las  can  grow  unstably  as  ona  nouns  back  in  tha  boon.  If  tha 
wuakar  condition  (37)  is  satisflad*  Instability  is  still  pradlctad*  but  of  an 


oscillatory  and  soaawhat  weaker  fora*  aoaauhat  dapandant  on  tha  aodal  and  tha 


▼alua  of  a.  Lao3  uses  a  -  0.7,  corresponding  to  a  Bennett  profile  truncated 
at  throe  to  four  Bennett  radii13.  Vo  can  write  tha  critical  value  of  Ip/l,ff 
for  instability*  froa  Bq.  (37)*  la  tha  fora  of  a  correction  to  (39), 


>  ffi2  ■  («  *  0.7*  X). 


Is 


Huaarlcal  evaluation  shows  that  H(a  •  0.7* X)  is  quite  does  to  unity*  varying 
only  froa  0.73  for  X  •  0  to  0.81  for  X  ♦  •*  so  (39)  is  a  reasonably  accurate 
instability  criterion  even  for  finite  a. 

Under  conditions  where  avalanche  is  uniaportant  even  at  tha  pinch  point* 
a.g.  for  a  typical  induction  llaec  baaa13*13  with  1^  <  10  ltd,  a#  ^  0.3  ca* 
in  air  at  density  close  to  or  above  aableat,  tha  Instability  condition  (41)  Is 
never  satisfied.  Ia  these  cnees*  Burp  sad  Leaps1  have  shown  that 
reaches  a  peak  value 


At' 


-r  '  ■ 
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I  /I  -  3  X  to  4  X 

p  «ff 


•t  the  pinch  point,  and  then  fells  off  to  •  fairly  constant  value 


I  /I  -  (1  ±  0.3)  X 

p  off 


(42a) 


(42b) 


a  few  cantina tars  further  back  in  (.  Equation  (42b)  is  incompatible  with 
condition  (41). 

If  avalanche  is  Important  near  the  pinch  point,  lp/ICf{  is  Increased. 
Hevertheleee  a  survey  performed  with  the  simulation  code  SMMO11,  which  does 
include  avalanche.  Indicates  that  over  a  very  wide  range  of  bean  parameters  in 


air  the  instability  condition  (41)  is  not  satisfied  over  a  long  enough  stretch 
of  beam  to  permit  effective  mode  growth.  [In  some  caaes  (41)  is  satisfied  for 
a  region  of  only  a  few  centimeters  about  the  pinch  point.]  For  exemple,  beams 
with  1^  -  10  kA  and  radius  £  0.5  cm  are  predicted  to  be  sausage-stable  in  air 
at  densities  above  50  torr,  and  beams  with  up  to  at  least  100  kA  are 
predicted  to  be  sausage-stable  in  air  at  ambient  density.  We  have  found, 
however,  that  the  instability  condition  (41)  is  satisfied  in  some  low-density 


regimes  where  avalanche  is  so  strong  near  the  pinch  point  that  X  ../X.  is 

til  D 

small,  e.g.  in  air  at  10  torr  X# ^ /X^  *  0.08  for  a  10  kA  beam  with  radius  0.5 
cm  and  current  risetime  0.33  nsec.  This  low  density  regime  around  10  torr  is 
of  Interest,  since  considerable  experimental  effort  has  been  devoted  to  it  in 
past20  and  present  experiments.** 

However  in  this  regime  avalanche  often  dominates  the  conductivity  physics 
even  well  behind  the  pinch  point.  Thus  the  conductivity  model  (5)  used  in  the 
present  theory  is  seriously  incomplete.  Inclusion  of  avalanche  eelf- 


eonelstently  in  the  model  of  the  been  body  probebly  would  be  stabilising. 


-ye- * *****  *4  ! 
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since  it  further  reduces  the  spetiel  seperetion  of  the  p Issue  current  frou  the 

been  current*  To  dete,  we  have  been  unable  to  treat  avalanche  self- 

consistently  in  an  analytic  theory,  but  we  are  presently  using  the 

axisyaaetric  beau  envelope  code  VIPER-0  to  siaulate  sausage  evolution  in  this 

regias.  In  Halted  studies  to  date  we  have  found  no  case  for  which  the 

21 

sausage  node  is  unstable  for  been  injection  into  neutral  gas. 

Recent  particle  ainulations8-10  have  shown  very  strong  axisyaaetric 
instabilities  in  sons  reglaes  where  (41)  predicts  stability.  We  present  a 
detailed  sluulatlon  analysis  of  these  instabilities  in  a  coapanion  paper1-1-. 

We  find  that  the  unstable  aode  is  the  (a  -  0,  n  «  2)  hollowing  aode,  not  the 
(a  •  0,  n  ■  1)  sausage  aode. 

For  beans  propagating  into  a  pre-lonlsed  channel,  can  be 

arbitrarily  large,  and  the  instability  conditions  (41)  or  even  (39)  can  be 
satisfied.  The  lialt  of  a  fixed  pre-lonlsed  channel,  l.e.  the  caae  in  which 
bean-induced  conductivity  augaantation  is  negligible  eoapared  to  the  pre¬ 
existing  conductivity,  has  been  considered  previously  by  Lee8.  In  our 
forasllsa,  this  is  the  Halt  A  ♦  0,  0  fixed,  A«(0)  fixed,  «(0)*»,  c 

reetrlcted  to  a  range  £  <  (  <  (  +&(,  where  u&c  «  l  and  o  (c)  has  the 

o  — '  —  o  o 

essentially  constant  value  <3.  Z  •  Since  <  is  proportional  to  A,  we  note 

no  o 

also  that  C  ♦  •  is  the  Halt.  The  dependence  of  any  perturbed  quantity  f 
o 

on  (,  froa  Bq.  (32),  then  reduces  to 
*  „  c  - 

f(c)  •  f  exp(-iw  ■ ■■■  —)  *  f  exp  [-ia(c  -  e  )],  (43) 

*o 

l.e.  exponential  growth  end  sinusoidal  oscillation  with  coaplex  frequency 


(45) 


Our  dispersion  relation  (34)  reduces  to  one  previously  derived  by  Lee 3, 

-a2  -  iaa  -  -  2  +  2  (%_  +  r-^L_r). 

T  Sff  Sff  1  55 


Our  condition  (39)  for  instability  to  occur  and  not  be  stabilised  by  any  value 
of  the  damping  coefficient  a  reduces  to 


I  /I  ~ 
p  eff 


>  3, 


(46) 


In  agreement  with  Ref.  5. 

In  conclusion,  then,  ve  have  determined  the  node  structure,  Eq.  (32),  for 

Instability  growth  in  the  bean  body,  and  have  calculated  a  dispersion 

relation,  Eq.  (34),  and  instability  conditions,  (39)  and  (41),  for  the  sausage 

■ode.  He  find  that  the  sausage  node  Instability  condition  is  not  satisfied 

for  beaas  injected  into  neutral  gas  under  conditions  that  satisfy  our 

assuaptlons.  However  the  instability  conditions  are  usually  satisfied  for 

beaas  injected  into  a  pre-ionlzed  gas  channel  with  the  saae  profile  as  the 

beaa  (and  with  4 wo  a  >  c),  because  of  the  substantial  current  neutralisation 
o  o  ~ 

under  those  conditions.  There  is  a  region  of  low  neutral  gas  density  where 
our  theory  would  predict  instability  but  where  the  theory  is  itself 
inapplicable  because  avalanche  is  strong  and  persistent.  Although  sausage 
stability  properties  are  not  well  understood  in  this  regias,  early  indications 
are  that  sausage  instability  does  not  occur. 

To  our  knowledge,  sausage  instability  has  not  been  observed 
experimentally  for  beams  injected  into  neutral  gas,  in  agreement  with  our 
conclusions. 
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